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We introduce a new mechanical model for producing tidally-driven strike-slip displacement along pre- 
existing faults on Europa, which we call shell tectonics. This model differs from previous models of 
strike-slip on icy satellites by incorporating a Coulomb failure criterion, approximating a viscoelastic rhe- 
ology, determining the slip direction based on the gradient of the tidal shear stress rather than its sign, 
and quantitatively determining the net offset over many orbits. This model allows us to predict the direc- 
tion of net displacement along faults and determine relative accumulation rate of displacement. To test 
the shell tectonics model, we generate global predictions of slip direction and compare them with the 
observed global pattern of strike-slip displacement on Europa in which left-lateral faults dominate far 
north of the equator, right-lateral faults dominate in the far south, and near-equatorial regions display 
a mixture of both types of faults. The shell tectonics model reproduces this global pattern. Incorporating 
a small obliquity into calculations of tidal stresses, which are used as inputs to the shell tectonics model, 
can also explain regional differences in strike-slip fault populations. We also discuss implications for fault 
azimuths, fault depth, and Europa’s tectonic history. 

© 201 1 Elsevier Inc. All rights reserved. 


1. Introduction 

Strike-slip offsets are common along Europan lineaments of 
varying length and type (Schenk and McKinnon, 1989; Tufts et al., 
1999; Hoppa et al., 2000; Figueredo and Greeley, 2000; Sarid et al., 
2002; Kattenhorn, 2002; Riley et al., 2006; Kattenhorn and Hurford, 
2009). A comprehensive survey of strike-slip faults in the Galileo 
regional mapping data set identified almost 200 strike-slip faults 
(Sarid et al., 2002). In Fig. 1, we show the azimuths and locations 
of the left lateral (red/dark gray) and right lateral (blue/light gray) 
faults that were identified. The survey revealed that, on a large-scale, 
the sense of motion along faults is not random. The distribution of 
faults is generally left lateral (LL) in the far north and right lateral 
(RL) in the far south, with a mixture of both fault types in between. 
The ratio of LL to RL faults in the mixed region decreases with lati- 
tude and trends differently between the leading and trailing hemi- 
spheres. Fault statistics from the survey (Sarid et al., 2002) are 
shown in more detail in Table 1. 

As Europa moves through its eccentric orbit, the magnitude and 
direction of its Jupiter-raised tidal bulges change. The result is dai- 
ly-varying deformation and stress that can drive tectonics. Tidal 
stress due to orbital eccentricity has been linked to the formation 
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of several types of tectonic features on Europa including strike-slip 
faults (Tufts et al., 1999; Hoppa et al., 2000; Kattenhorn, 2002; Sarid 
et al., 2002; Rhoden et al, 2011). The motion and eruption timing 
along the tiger stripe fractures on Enceladus (Hurford et al., 2007; 
Nimmo et al., 2007; Smith-Konter and Pappalardo, 2008) have also 
been attributed to tides caused by Enceladus’ eccentricity and close 
proximity to Saturn. 

Additional orbital and rotational characteristics of these satel- 
lites can also influence the pattern of stress change throughout 
each orbit. Gravitational interactions between Jupiter’s large satel- 
lites cause Europa’s obliquity to be at least 0.1°, and possibly larger 
if the satellite has a subsurface ocean (Bills et al., 2009). Evidence of 
~1° of obliquity has been found within Europa’s geologic record, 
preserved in the shapes of arcuate features called cycloids (Hurford 
et al., 2009; Rhoden et al., 2010). Rhoden et al. (2010) found that 
good fits to cycloids were achieved only when longitude transla- 
tion was assumed, supporting the idea that Europa’s shell rotates 
non-synchronously (Rhoden et al., 2010). However, cycloid fits 
did not indicate that stress from non-synchronous rotation (NSR) 
influences their formation. Non-synchronous rotation of Europa’s 
ice shell has been proposed on theoretical grounds (Greenberg 
and Weidenschilling, 1984; Ojakangas and Stevenson, 1989) and 
is supported by geologic evidence (see review in Bills et al. 
(2009)). However, if the rotation is slow, the accumulation of stress 
from non-synchronous rotation may be limited (Greenberg et al., 
1998; Goldreich and Mitchell, 2010). 
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Observed Strike Slip Sense 



Longitude 

Fig. 1. Observed faults from the survey by Sarid et al. (2002) shown at their 
measured azimuths. Left-lateral faults are blue (dark gray), and right-lateral faults 
are red (light gray). The faults are binned in latitude and longitude. At longitude 
90°W and latitudes 0° and 15°S, pervasive chaos formation inhibited fault 
identification. Other empty circles represent regions that were not mapped due 
to limited high-resolution imagery. (For interpretation of the references to color in 
this figure legend, the reader is referred to the web version of this article.) 


Table 1 

Strike-slip statistics from Sarid et al. (2002). 



Trailing hemisphere 

Leading hemisphere 

LL 

RL 

% LL 

LL 

RL 

% LL 

60N+ 

18 

0 

100 

- 

- 

- 

40-60N 

8 

6 

57 

2 

0 

100 

20-40N 

17 

13 

57 

4 

1 

80 

0-20N 

1 

4 

20 

3 

2 

60 

20S a -0 

7 

47 

13 

- 

- 

- 

40-20S 

0 

8 

0 

3 

5 

38 

60-40S 

0 

11 

0 

8 

16 

33 

60S- 

0 

7 

0 

0 

1 

0 


a The southernmost LL fault in the trailing hemisphere is located at 21.4°S; we 
included it in the 0-20°S bin. 


The tidal stress due to orbital eccentricity and obliquity can be 
decomposed into normal and shear components along a fault of 
given azimuth and location (Section 3.1). Depending on the fault 
response model adopted, the changes in stress throughout an orbit 
can be analyzed to determine the likely net offset direction along 
the fault. We present a new mechanical model for the formation 
of these tidally-driven strike-slip faults (Sections 2.2 and 3.2). 
Our “shell tectonics” model combines tidal stress and an elastic re- 
bound seismic model for slip and stress release on faults, in which 
accumulated shear stress on the fault is released when the fault 
slips. Our model successfully reproduces the observed global pat- 
tern of strike-slip faults on Europa. When obliquity is also incorpo- 
rated, we can generate variations in the global pattern that are 
consistent with the observed differences in the leading and trailing 
hemisphere strike-slip fault populations (Sections 4 and 5.1). In 
contrast, we find that a model neglecting the release of stress 
due to slip is less effective at matching the global pattern of 
strike-slip (Sections 4 and 5.1). We also discuss implications of 
these results on relative offset accumulation rates along faults 
and the depth to which strike-slip motion may occur (Sections 
5.2 and 5.3). In addition to advancing our understanding of 
tidally-driven tectonics on Europa, the shell tectonics model may 
also be relevant for Enceladus. Strike-slip motion and subsequent 
heat generation along the Tiger Stripe fractures has been linked 


to tides (Nimmo et al., 2007; Smith-Konter and Pappalardo, 
2008), and tidally-driven opening and closing of the fractures has 
been proposed as a mechanism for controlling plume eruptions 
(Hurford et al., 2007). 

2. Models of tidally-controlled strike-slip 

2.1. Previous models 

As Europa moves through its eccentric orbit, a pre-existing fault 
will experience diurnally-varying normal and shear stresses that 
may drive tectonics. Hoppa et al. (1999, 2000) introduced the tidal 
walking model, which assumes that faults slip freely while in ten- 
sion, but that slip is prohibited during compression. Once failure 
initiates, slip on the fault accumulates until the magnitude of the 
tidal shear stress reaches a maximum. Then, because the ice shell 
is assumed to behave elastically over short (i.e. diurnal) timescales, 
slip on the fault is reduced back toward zero as the shear stress de- 
creases. Therefore, if the shear stress crosses zero during the ten- 
sion phase, it is assumed that any accumulated slip up to that 
point has been recovered, and there is no net offset on the fault. 
When the normal stress becomes compressive, the fault is 
clamped; an offset on the fault reflects the amount of accumulated 
slip since the last time the shear stress crossed zero. In order to 
maintain this offset over subsequent orbits, elastic stresses in the 
shell must relax while the fault is clamped so that, when the nor- 
mal stress again becomes tensile, the fault will not slip back to its 
initial configuration. 

Hoppa et al. (1999, 2000) used the tidal walking model to create 
global predictions of slip direction, which generally match the ob- 
served global pattern on Europa. Rhoden et al. (2011) made a cor- 
rection to the predictions, but this did not affect the resulting 
global pattern of only left lateral (LL) faults poleward of 35°N, only 
right lateral (RL) faults poleward of 35°S, and between these re- 
gions, either right or left lateral faults with the slip-direction 
depending on the longitude and the azimuth of the crack. Initially, 
observed differences in the strike-slip statistics of the leading and 
trailing hemispheres were attributed to non-synchronous rotation 
and true polar wander (Sarid et al., 2002). However, Rhoden et al. 
(2011) investigated the effects of obliquity on the predictions and 
concluded that an obliquity of ~1° could lead to offsets of the re- 
gions of mixed right and left lateral faults that are more consistent 
with the observations than polar wander. 

Its ability to reproduce the observed global pattern of strike-slip 
displacement is a strong indication that the tidal walking model 
captures a fundamental attribute of the fault slip process on Euro- 
pa. However, the mechanics of the model are simplified. Neither 
stress release due to slip nor stress relaxation with time are explic- 
itly included in the tidal walking model. In fact, predictions of slip 
direction are based solely on the sign of the shear stress when the 
normal stress becomes compressive (see Rhoden et al., 2011). Fur- 
thermore, failure is only allowed when the fault experiences ten- 
sile normal stress so the model is only applicable in the upper 
tens of meters of Europa’s surface. Deeper within the shell, the 
compressive overburden stress would exceed the tidal stress and 
clamp the fault throughout the entire orbit. 

A model for strike-slip formation due to tidal stress has also 
been proposed for Enceladus (Smith-Konter and Pappalardo, 
2008) and should be equally applicable to Europa. In this model, 
which we will refer to as the SKP model, the fault can slip when- 
ever the tidal shear stress exceeds the Coulomb failure criterion. 
Slip on the fault accumulates via a large event when the shear 
stress exceeds the Coulomb failure threshold and subsequent 
small, creeping events that continue as long as the tidal shear 
stress remains above the failure threshold. Summing the accumu- 
lated slip over one orbit yields a prediction of slip direction. 
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Comparison between model predictions and the observed tiger 
stripe fractures on Enceladus is inconclusive (Smith-Konter and 
Pappalardo, 2008), and the model has not yet been tested against 
the global pattern of strike-slip faults on Europa. 

A major difference between the tidal walking model and the SKP 
model is in the determination of slip direction throughout an orbit. 
In the tidal walking model, it is assumed that the shell, being elastic 
over short timescales, will act to restore the fault to its initial config- 
uration. Hence, slip depends on the change in shear stress with time. 
In the SKP model, fault slip is determined by the value of the tidal 
shear stress regardless of whether it is increasing or decreasing. 
One way to visualize the difference is to think of the SKP model as 
describing the behavior along a boundary between two uncon- 
nected plates. Whereas, in the tidal walking model, the plates are 
connected in the far field by a spring. It is unclear a priori , which 
model best represents the behavior of faults on Europa and 
Enceladus. 

To further investigate the mechanism of strike-slip fault forma- 
tion on Europa, we take a dual approach. First, we develop a model 
based on the assumption of elastic rebound over short timescales, 
in the spirit of the tidal walking model, but that also includes a more 
quantitative and physical treatment of stress accumulation and slip 
on a fault. We also create a simplified version of the SKP model, 
adopting the rules for stress accumulation and slip outlined by 
Smith-Konter and Pappalardo (2008) and developed with their guid- 
ance (Smith-Konter, personal communication, 2010). We then test 
the underlying fault slip assumptions of each model by creating pre- 
dictions of fault slip direction and comparing the predictions against 
the observations in the strike-slip survey by Sarid et al. (2002). 

2.2. Shell tectonics 

In the shell tectonics model, slip occurs when the accumulated 
shear stress on the fault exceeds the Coulomb failure threshold. 
When a fault is in tension, we set the failure criterion to zero (no 
cohesion). This results in continuous slip while the fault is in tension 
and large discrete events in compression. As with the SKP model for 
Enceladus (Smith-Konter and Pappalardo, 2008), allowing slip while 
the fault is under compression allows the shell tectonics model to 
operate deeper in the ice shell than the tidal walking model. 

When the failure threshold is exceeded, the amount of resulting 
slip is assumed to be proportional to the accumulated shear stress 
on the fault at that time. We explicitly include the release of accu- 
mulated shear stress due to slip by resetting the shear stress on the 


fault to zero after every slip event. Hence, the tidal shear stress has 
to increase in magnitude for the fault to continue slipping in the 
same direction. If the shear stress instead decreases, the fault will 
begin to slip back toward its initial configuration. Thus, release of 
accumulated stress accounts for the elastic response of the shell. 
This is in contrast to the SKP model in which the fault continues 
to creep after the failure threshold is reached, and the direction 
of creep is always in the same direction as the initial slip event. 
To account for the long-term viscous behavior of Europa’s ice shell, 
we also allow a small amount stress of relaxation when on the 
fault, which we find to be critical to producing net offsets. To make 
a prediction of slip direction, we determine the long-term trend of 
motion on the fault by calculating the difference in the net offset 
over many orbits. Additional details of these calculations are given 
in Section 3.2. 

Use of accumulated shear stress rather than tidal shear stress 
when determining failure is a keystone of the shell tectonics mod- 
el. Fig. 2 is a schematic showing the tidal deformation of an uns- 
lipped region and a fault’s response to that deformation at 
certain points in the orbit. The direction of the most compressive 
principal stress is also shown. Throughout the schematic, there is 
a regional left lateral (LL) tidal shear stress, which is increasing 
in magnitude at A and B but decreasing thereafter. We also assume 
that the regional tidal normal stress has just entered a period of 
compression. Previously, the fault was free to slip in response to 
the left lateral tidal deformation because it was in tension. 

At A, the fault has a LL offset and zero accumulated shear stress. 
At B, the tidal deformation in the region has become increasingly 
left lateral. Based solely on the regional deformation, the fault 
would accumulate another LL offset, even larger than at A. How- 
ever, if we only consider the accumulated stress on the fault, we 
get a much different answer. Because, the fault has already re- 
sponded to the deformation at A, it accumulates a small amount 
of additional shear stress due to the incremental increase in defor- 
mation, but it is not enough to reach the failure threshold. At C, the 
tidal shear stress has decreased and is now identical to A. Thus, the 
accumulated shear stress on the fault returns to zero. At D, the re- 
gional deformation has decreased further. If the fault had never 
slipped, we would expect the deformation to induce left-lateral 
slip on the fault as long as the tidal stress exceeded the failure cri- 
terion. However, the fault has already slipped to accommodate the 
deformation at C, which is now more LL slip than would be gener- 
ated at D. Hence, at D, the fault actually accumulates right lateral 
shear stress. At E, the deformation has decreased so much that 


LL shear stress increasing 





LL shear stress decreasing 

C D E 


[S 1 0 0 





Fig. 2. A schematic of tidal deformation and fault response with changing tidal shear stress. The direction of the most compressive principal stress is also shown. In A and B, 
the shear stress is increasing generating left lateral slip on the fault and reducing the accumulated shear stress on the fault. In C-E, the shear stress is decreasing. Because 
stress has already been relieved by slip, the fault begins to accumulate right lateral shear stress and ultimately slips in that direction, thereby reducing the net offset on the 
fault. 
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the accumulated shear stress on the fault generates right lateral 
slip even though the regional tidal shear stress is still left lateral. 
The shell tectonics model would thus predict left lateral slip at A 
and right lateral slip at E. In the SKP model, the fault would under- 
go a large LL slip at A, continue creeping in a LL sense until C, and 
not slip again in this time period because the tidal shear stress is 
lower than the failure threshold. In the tidal walking model, no slip 
would have occurred during this entire time period because the 
fault was in compression. 

3. Generating global predictions 

3.2. Calculating tidal stress 

Our formation model is based on the hypothesis that tidal 
deformation is the driver of strike-slip displacement on Europa. 
The daily changes in tidal stress thus control the ability of pre- 
existing faults to slip and the direction of net shear displacement 
along the faults. We neglect stress from non-synchronous rotation 
or overburden pressure (see also, Section 5.3). We assume that 
Europa has long-since adopted its primary tidal shape, which is 
based on the average distance and angle between Europa and Jupi- 
ter. It is deviations away from that shape that generate stress. We 
use the equations for stress in a thin, elastic shell (e.g., Melosh, 
1977, 1980): 

g 3 = C(5 + 3 cos 2 Sp) (la) 

Go, = -C( 1 - 9 cos 2 Sp) (lb) 

where C = 3h 2 Mp{\ + v)l8npa 3 (S + v) and h 2 is Europa’s tidal Love 
number, M is Jupiter’s mass, p is Europa’s shear modulus, v is the 
Poisson’s ratio of Europa, p is Europa’s average bulk density, a is 
the average distance between Jupiter and Europa, and S P is the 
angular distance from a point on Europa’s surface to its primary ti- 
dal bulge. The a§ stress is directed radially from the tidal bulges and 
the G a stress is perpendicular to a d . Because all the stress calcula- 
tions contain a factor of C, its value does not influence our results. 

The tidal stress equations can be modified to account for Euro- 
pa’s eccentricity and obliquity (Hurford et al., 2009; Rhoden et al., 
2010, 2011). Eccentricity causes the tidal bulges to librate in longi- 
tude while obliquity predominantly causes a latitudinal libration. 
Spherical trigonometry is used to calculate the time-varying loca- 
tion of the tidal bulge, which now depends on the spin pole direc- 
tion and the true anomaly. Because the locations of the tidal bulges 
are changing, the angular distance to the bulge, < 5 , also changes. The 
equations then become: 


g$ = C(1 — ecosn) 3 (5 + 3cos2<S) 

(2a) 

G a = -C( 1 - ecosn) _3 (l - 9cos2<S) 

(2b) 

and 


Bulge colatitude = 7i/2 - 8sin(n + cp) 

(2c) 

Bulge longitude = -2esinn 

(2d) 


where e is Europa’s eccentricity, s is Europa’s obliquity, (p is the spin 
pole direction (SPD), and n is the true anomaly. When Europa is at 
pericenter, if the spin pole is pointing toward Jupiter, the SPD is de- 
fined as 90°; SPD increases clockwise. There is degeneracy between 
spin pole direction and longitude such that the stress field is iden- 
tical when both are modulated by 180°. 

Stress from the diurnal tide is calculated by subtracting the pri- 
mary tidal stress (Eqs. (la) and (lb)) from the stress due to both 
eccentricity and obliquity (Eqs. (2a) and (2b)) once the stresses 
have been rotated to a common coordinate system. The resulting 
diurnal stresses in the new coordinate system are in the direc- 
tion of the north pole, and the perpendicular stress, g a * There may 


now also be a shear stress, o d0L * These diurnal tidal stresses can 
then be decomposed into normal and shear components relative 
to a fault’s orientation, where £ is the azimuth of the crack mea- 
sured clockwise from north. 

g = 0.5((J^ + <7 a *) + 0.5(<7«5* - (7a*) cos(2£) + a**. sin(2Q (3a) 
t = -0.5((7*, - g a *) sin(2Q + g 5oc * cos( 2Q (3b) 

To determine the slip direction along a fault, we first specify the 
latitude and longitude of the fault, the fault azimuth, and the pre- 
scribed amount of obliquity and spin pole direction. Using these 
parameters, we calculate the stress throughout an orbit. 

The equations for tidal stress shown here assume an elastic 
rather than viscoelastic shell. Wahr et al. (2009) present equations 
for tidal stress in a viscoelastic shell. However, the Wahr et al. 
(2009) model does not include obliquity, an important parameter 
in our work. In addition, the predictions made with the tidal walk- 
ing model relied on the formulation we present here (Hoppa et al., 
1999, 2000; Rhoden et al., 2011), so it is useful for consistency to 
test shell tectonics against the same stresses. Smith-Konter and 
Pappalardo (2008) did use the Wahr et al. (2009) model, so it is 
important to consider how the two approaches differ in their cal- 
culations of tidal stress. 

A major result of Wahr et al. (2009) is that stresses caused by 
non-synchronous rotation in a viscoelastic shell are phase-shifted 
from those in an elastic shell. For diurnal stresses, the amplitudes 
can differ by 0.1-3% between the elastic and viscoelastic cases, 
depending on the viscosity assumed for the ice shell. Because we 
neglect the effects of NSR stress, the differences between the diur- 
nal stresses we calculate here and those used by Smith-Konter and 
Pappalardo (2008) should be small - indeed, no larger than 3%. Our 
main goal is to test whether observations are best explained with 
the assumption that - quite generally - slip on a fault is reduced 
when the shear stress decreases (shell tectonics) or that slip con- 
tinues to increase even if the shear stress decreases (SKP). The ex- 
act amplitude of the tidal stress will not significantly affect this 
determination. 

3.2. Applying shell tectonics 

We use the shell tectonics model to fit the global strike-slip pat- 
tern on Europa by generating predictions on a global scale. We 
determine net offsets for faults with azimuths from 0° to 180° in 
1° increments, at longitudes from 0° to 360° every 30°, and from 
latitude 75° to -75° in increments of 15°. We apply shell tectonics 
to tidal stress fields with and without 1° of obliquity and test two 
end-member values of the friction coefficient: 0.6 and 0.2. The 
higher friction value is appropriate for rock and may be applicable 
to very cold ice (Beeman et al., 1988); the lower value was used by 
Smith-Konter and Pappalardo (2008). We use it here for 
consistency. 

To calculate the accumulated stress on the fault, we solve the 
following differential equation in which a tida i is the stress tensor 
and stress relaxes on an e-folding timescale of I/77: 

dO a cc/dt = (dffthu/dt) - Y}<5 acc (4) 

The relaxation rate is a depth-averaged value that accounts for the 
decrease in ice viscosity with depth. Our solution is restricted to 
values of rj for which pt 0 r bit < 1. Without relaxation, we find that 
no net slip accumulates over time. 

There is no release of accumulated normal stress in this model. 
Rather, we assume that the accumulated normal stress on the fault 
reflects the long-term behavior of daily-varying stresses. We ana- 
lytically solve Eq. (4) to determine the initial values of the accumu- 
lated normal stress components such that the transient part of the 
solution is minimized for each component. We also advance in the 
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orbit to a time when the fault is in tension so the initial value of the 
accumulated shear stress can be set to zero. At each subsequent 
time step, we first add Aa t idai and calculate the normal and shear 
components of the accumulated stress (using Eqs. (3a) and (3b)). If 
the accumulated normal stress on the fault is positive (in tension), 
it will slip for any non-zero value of accumulated shear stress. In 
compression (a acc < 0), the accumulated shear stress must exceed 
the Coulomb failure threshold. 

^acc — —f^acc (5) 

where /is the friction coefficient. If slip occurs, the cumulative slip 
is increased by the amount of accumulated stress on the fault and 
the accumulated shear stress is set to zero. Finally, we reduce the 
deviatoric component of any remaining accumulated stress, 
Gacc. according to Eq. (4) to account for the long-term viscoelastic 
behavior of the shell. We compute the accumulated stress and 
cumulative slip at 8500 time steps per orbit and run each simula- 
tion over 1000 orbits to ensure that the solution has converged. 

To predict the slip direction on the fault over time, we calculate 
the difference in slip at the same point in two successive orbits to 
determine the net offset. We then average the net offset over 1000 
orbits to determine the long-term behavior on the fault. To reduce 
numerical noise, we repeat the simulation 20 times with a random 
number (between 0 and 1 ) added to the longitude, and average the 
net offset values over the 20 simulations. If the average net offset is 
positive, we predict left-lateral displacement on the fault; a nega- 
tive average net offset yields a right-lateral prediction. If the net 
offset averages to zero, we predict no net slip on the fault. 

3.3. Applying an SKP-type model 

Because the observed global pattern is such a useful metric for 
testing tidal-tectonic models, we also generate predictions using a 
model based on Smith-Konter and Pappalardo (2008). In that model, 
there is a large slip event when the fault first reaches the failure cri- 
terion. Then, small slips occur as long as the tidal stress on the fault 
exceeds the Coulomb failure threshold. The amount of slip from the 
large events is proportional to the tidal shear stress on the fault at 
the time of failure. For the creeping events, Smith-Konter and Pap- 
palardo (2008) assume a constant, time-averaged, strain rate to 
which the amount of accumulated slip is proportional. Unlike tidal 
walking or shell tectonics, the direction of slip in the SKP model is 
always the same as the tidal shear stress direction. The actual model 
presented in Smith-Konter and Pappalardo (2008) is more sophisti- 
cated than what we are using here. However, our interest is in test- 
ing the determination of when and how fault slip will occur in the 
SKP model versus the shell tectonics model. To accomplish this, 
we have adopted a set of rules for stress accumulation and fault slip 
that capture these fundamental assumptions. 

To compute a prediction of direction using our SKP-type model, 
we keep a running sum of LL (positive) and RL (negative) slip 
throughout an orbit. When the tidal stress on a fault first reaches 
the failure criterion for one slip direction, there is a large slip event 
in which all of the tidal stress on the fault is relieved. Hence, the 
amount of slip is proportional to the amount of shear stress. Sub- 
sequently, we calculate the incremental change in stress at each 
time step within the slip window; the rate of change of the stress 
is proportional to the slip rate and can thus be used to determine 
additional slip. Eventually, the shear stress becomes lower than 
the failure threshold and no slip occurs on the fault. Later, the fault 
may slip in the opposite direction reducing the total offset on the 
fault. After one orbit, positive cumulative slip would mean a LL pre- 
diction for the fault, negative would be RL, and zero net slip implies 
no offset on the fault. In practice, we find that the failure stress in 
each slip direction is similar so the prediction is most dependent 
on the length of the slip window in each direction. Only in cases 



Fig. 3. Shell tectonics predictions of slip direction with zero obliquity. Within each 
circle, black regions indicate crack azimuths along which we predict left lateral 
displacement; light gray represents right lateral fault azimuths. We predict no net 
slip at the boundaries between left and right regions. The last column shows the 
predictions summed over all longitudes, in which dark gray represents azimuths 
that could have right or left lateral displacement depending on their longitude at 
the time the displacement occurred. The pattern of slip direction generated using 
the shell tectonics model fits with the global observed pattern, in which left-lateral 
faults dominate in the far north, right lateral faults dominate in the far south, and 
there is a mixture of right and left lateral faults in between. 



Fig. 4. Shell tectonics predictions of slip direction with 1.0° obliquity at a spin pole 
direction of 220°. Black represents fault azimuths that will exhibit left lateral 
displacement, and light gray represents right lateral faults. The last column shows 
the predictions at all longitudes combined, in which dark gray indicates that either 
type of fault can form over the span of longitudes at which we calculated 
predictions. Obliquity breaks the longitudinal symmetry in the pattern of predicted 
slip direction, which is consistent with the observed differences in the observed 
fault populations of the leading and trailing hemispheres. 

in which the fault spends equal time in LL and RL slip does the 
amount of initial slip play a role. This result is consistent with 
the results of Smith-Konter and Pappalardo (2008). 

4. Results 

The shell tectonics global predictions accounting for only eccen- 
tricity are shown in Fig. 3. The predictions using tidal stresses that 
include a 1° obliquity and spin pole direction of 220° are shown in 
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Fig. 4. Both sets of predictions use the lower friction value of 0.2; 
predictions using the higher friction value are shown in the Elec- 
tronic supplement. In these plots, black represents fault azimuths 
that undergo left lateral slip, and light gray represents right lateral 
slip. The net slip goes to zero at transitions between left and right 
lateral regions. In the last column of each plot, we show a combi- 
nation of the predictions made at each latitude over all longitudes. 
Here, dark gray indicates that both right and left lateral faults form 
within the span of longitudes tested for a given crack azimuth. To 
generate the predictions shown in Figs. 3 and 4, we use a relaxation 
rate of 1 x 10 -10 s -1 , which is comparable to the Maxwell time of 
ice at ~130 K. The relaxation rate will determine the magnitude of 
the accumulated offset per orbit. However, our predictions of slip 
direction and relative offset accumulation are insensitive to the 
relaxation rate over several orders of magnitude (see online Sup- 
plement for an example). 

Without obliquity (Fig. 3), the shell tectonics model predicts re- 
gions with only left lateral faults at latitudes 60° and 75°, and right 
lateral faults at -60° and -75. Latitudes 45° and -45° are nearly 
exclusive, but there are a few longitudes with very narrow azimuth 
ranges that predict the opposite slip direction. At longitude 90°W 
and 270°W, on the other hand, the regions of exclusivity extend to- 
ward the equator to include ±30°. At all other latitudes, there is a 
mixture of RL and LL faults predicted at every longitude. Using 
the higher friction value (0.6) produces very similar global predic- 
tions for this case, as shown in the Electronic supplement. 

With obliquity, a distinct hemispheric pattern emerges, as 
shown in Fig. 4 (for (p = 220°). In the trailing hemisphere, around 
longitude 210°W, exclusively LL regions are predicted farther north 
than without obliquity while exclusively RL regions are predicted 
closer to the equator. In the leading hemisphere (Ion 90°W), the 
opposite occurs: LL-only regions begin farther south, and mixed re- 
gions extend into lower latitudes than without obliquity. With the 
higher friction value (see SOM), there are similar trends, but fewer 
locations are predicted to have exclusively LL or RL faults for all 
azimuths. The spin pole direction determines the longitudes at 
which the asymmetries in predictions would occur. Qualitatively, 
the pattern of offsets produced by obliquity appears consistent 
with the hemispheric differences in the observed strike-slip 
populations. 



Fig. 5. Model predictions of slip direction with zero obliquity using the SKP-type 
model. Black represents fault azimuths that will exhibit left lateral displacement, 
and light gray represents right lateral faults. The last column shows the predictions 
at all longitudes combined, in which dark gray indicates that either type of fault can 
form over the span of longitudes at which we calculated predictions. This model 
predicts roughly equal numbers of right and left lateral faults at any location given a 
random distribution of fault azimuths. 


The predictions using the SKP-type model (without obliquity) 
are shown in Fig. 5. In all regions, both left and right lateral faults 
are predicted to form with the slip direction depending on fault 
azimuth. Given a population of faults with random azimuths, this 
model should produce equal numbers of LL and RL faults at all lat- 
itudes and longitudes on Europa. 

5. Discussion 

5.1. Matching the slip directions of observed faults 

To quantify the goodness of fit between the shell tectonics pre- 
dictions and the observed faults, we made predictions of slip direc- 
tion at the actual latitude, longitude, and azimuth of each observed 
fault (as reported in Rhoden et al., 201 1 ). As shown in Figs. 3 and 4, 
the hemispheric differences in the observed fault population are 
better matched when obliquity is included in the tidal stress calcu- 
lations. By testing predictions made with different spin pole direc- 
tions against the individual observed faults, we find that an SPD of 
220° produces a pattern that is most consistent with the observed 
leading-trailing asymmetry. Using high friction, the slip directions 
of 70% of faults in the survey by Sarid et al. (2002) are correctly 
predicted at their current longitudes. Using the lower friction va- 
lue, our fit improves to 75%. Rhoden et al. (2011) tested the tidal 
walking model predictions against the same survey of faults. They 
found that an obliquity of 1° and SPD of 270° could account for the 
slip directions of 69% of faults - their best fit without appealing to 
longitude migration. The underlying assumptions concerning fault 
slip are similar between the shell tectonics model and the tidal 
walking model. However, Rhoden et al. (2011) only tested every 
90° of SPD whereas we tested at every 10°; that may account for 
the difference in optimal spin pole direction. The SKP-type model 
correctly predicts the slip directions of only 46% and 50% of the ob- 
served faults with eccentricity-only and with 1° of obliquity 
{cp = 220°), respectively. Because the SKP-type model does not pro- 
duce a hemispheric asymmetry, the value of the spin pole direction 
does not significantly influence the goodness of fit to the 
observations. 

In Fig. 6, we show the regions included in the strike-slip survey 
by Sarid et al. (2002), with different colors representing the per- 
centage of fault slip directions accurately predicted with the shell 
tectonics and SKP-type models. The total number of faults in each 
region is also listed on shell tectonics model graphic. The numbers 
are the same for the SKP-type model. The shell tectonics model 
performs well in most regions but has trouble in some areas of 
mixed LL and RL faults. The SKP-type model predicts mixed faults 
in all regions, so it performs best in mixed regions and worst in re- 
gions where one slip direction is observed more than the other. 

5.2. Matching the azimuth distribution of observed faults 

Using shell tectonics, we calculate the net fault offsets that are 
predicted to accumulate during each orbit. The absolute magni- 
tudes of these offsets depend on poorly-constrained model param- 
eters, so we instead compare the relative accumulation rates along 
faults. We find that, within a given region, faults with certain azi- 
muths are predicted to develop larger offsets per orbit than others. 
Faults at these azimuths should be easier to identify as having 
strike-slip displacement. The shell tectonics model thus predicts 
differences in observed fault azimuth distributions in different re- 
gions based on the relative accumulation rates. Some regions with- 
in the strike-slip survey do imply a preferred direction, such as the 
left lateral faults at 75°N that trend mostly east-west (Fig. 1 ). How- 
ever, fault at most azimuths are observed in the equatorial regions 
at longitude 210°. 
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Fig. 6. The regional accuracy of the predictions of fault slip direction using the shell tectonics model (left) and SKP-type model (right) with an obliquity of 1.0° and a spin pole 
direction of 220°. The shading represents the percentage of faults with slip directions that were accurately predicted with each model. The number of faults in each region is 
shown in the shell tectonics model graphic; the numbers are the same for the SKP-type model. Regions that are not shaded contain no observed faults. 
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Fig. 7. Relative rates of offset accumulation for a given region using the shell 
tectonics model (obliquity 1°, SPD 220°) are shown in an annulus surrounding the 
observations from Fig. 1. The gradient of rates goes from fast accumulation of left 
lateral slip (bright red) to fast accumulation of right lateral slip (bright blue). White 
can represent slip in either direction but at a reduced pace compared to other 
azimuths. The differences in accumulation rate may give rise to an observation bias 
within the fault population. Faults that accumulate slip faster than neighboring 
faults should be easier to identify. (For interpretation of the references to color in 
this figure legend, the reader is referred to the web version of this article.) 
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We apply a statistical test to determine the extent to which the 
shell tectonics model can account for the azimuth distribution of 
the observed strike-slip fault population. In Fig. 7, we show the 
observed fault azimuths from Fig. 1 embedded within an outer 


annulus that displays the relative offset accumulation rate in the 
region from the shell tectonics model. Red represents the fastest 
(or largest) accumulation of LL slip, blue shows the fastest (or larg- 
est) accumulation of RL slip, and white regions may produce either 
LL or RL faults (as shown in Fig. 4) but at a reduced rate compared 
to other azimuths. These rates were calculated using the parameter 
values that provide the best fit to the global pattern (e = l°, 

(p = 220°,/= 0.2). 

Fig. 7 shows that in some regions there should be preferred ori- 
entations at which offsets grow more rapidly. We can calculate the 
probability of producing the observed distribution of fault azi- 
muths with the model shown in Fig. 7, as follows. For a given re- 
gion (one circle in Fig. 7), we can construct a curve of the offset 
accumulation rate with azimuth. Normalizing the curve converts 
the rates to probabilities. Fig. 8 shows the probability curve for lat- 
itude 75°N and longitude 210°W as an example. In this region, 
faults of all azimuths are predicted to accumulate LL offsets; the 
peak in probability from ~80° to ~135° yields a prediction that 
more faults will be observed at these azimuths. Following this ap- 
proach, it is straightforward to calculate the probability of forming 
each fault in the observed strike-slip population based on its bin- 
ned latitude and longitude and its azimuth, rounded to the nearest 
degree. To compute a probability for all the faults in the survey, the 
probabilities for each fault are multiplied. The SKP-type model also 
yields net offsets along faults, so we use the same approach to test 
the probability of that model producing the observed fault 
population. 

Neither the shell tectonics model nor the SKP-type model cor- 
rectly predicts the slip directions of all of the observed faults. A 
fault with a slip direction that is not accurately predicted has zero 
probability of forming in the model. The probability of the model, 
given the observations, is proportional to the product of the indi- 
vidual fault probabilities, so we must either exclude the faults with 
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Azimuth 

Fig. 8. The probability versus azimuth curve for 75°N, 210°W based on the relative offset accumulation rates predicted with the shell tectonics model for an obliquity of 1.0° 
and a spin pole direction of 220°. 



Fig. 9. The probabilities of forming the observed fault population using a two-component model for the shell tectonics (circles) and SKP (squares) models. In the two- 
component test, we assume that the slip directions and azimuths at which some offsets form are determined by the shell tectonics or SKP-type model while others are formed 
with random slip directions and azimuths. In the highest probability model, half the observed faults are formed via shell tectonics and half are formed randomly. The highest 
probability model that uses the SKP-type model requires that 75% of the observed faults formed at random orientations. 


slip directions that are not accurately predicted or include a popu- 
lation of faults that formed with random orientations. For a fault 
had an equal likelihood of forming at any azimuth and either slip 
direction, the probability of generating that fault would be (1/ 
360). Generating a specific set of randomly-selected faults would 
have a probability of (l/360) n , where n is the total number of 
faults. 

We first test the probability of forming the observed population 
with a two-component model, in which some faults formed via the 
shell tectonics or SKP-type model and some formed at random ori- 
entations through some other process. The probability of forming 
each observed fault, P obs , is: 

P 0bs =/r/360 + (1 -fr) * P mod el (6) 

where f T is the fraction of faults that formed at random orientations 
and P m odei is the probability of forming the fault with the model 
(either shell tectonics or SKP-type) using the offset accumulation 
rates described above. As shown in Fig. 9, the two-component mod- 
el with the highest probability of forming the observed fault popu- 
lation is one in which half the faults were created via shell tectonics 
and half were created at random azimuths. The SKP-type model 


performs best when 75% of the faults are part of the random set; 
the probability of this model creating the observed fault population 
is lower than the best overall model by a factor of 10 20 . 

With both models, as the fraction of randomly-oriented faults 
approaches zero, the probability of either model forming the ob- 
served population decreases because the faults with the wrong 
predicted slip direction can no longer be accounted for. The popu- 
lation of randomly-oriented faults may reflect a population of 
faults that have migrated in longitude since their formation. In 
addition, a region may contain faults with different ages. An older 
fault could have built up enough offset over time to appear more 
prominent than a young fault even if the young fault could develop 
a net offset more quickly. A history of offset development and lon- 
gitude migration could be pieced together with enough sequence 
information. The fault history would then reflect the rate of non- 
synchronous rotation of the surface and the rate of growth and 
destruction of strike-slip offsets. Unfortunately, limited imagery 
makes sequencing a difficult task. Many areas simply do not have 
enough unambiguous cross-cutting relationships to fully constrain 
the formation sequence (e.g. Sarid et al., 2006). More complete 
sequencing should be possible if additional imagery could be 
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Fig. 10. The regional ratios of the probability of forming the observed faults using 
the shell tectonics model to that of a model that assumes random orientations and 
slip directions. This analysis includes only those faults with slip directions that were 
accurately predicted with the shell tectonics model; the number displayed within 
each region is the number of correctly predicted faults. Regions in shades of green 
represent regions in which the shell tectonics model outperforms the random 
model. Regions in yellow contain only one fault, so their usefulness is limited. The 
area shown in red is the one region that has more than one fault and is better 
explained by the random model than the shell tectonics model. (For interpretation 
of the references to color in this figure legend, the reader is referred to the web 
version of this article.) 



obtained with at least the ~250 m/pix resolution of the regional 
mapping imagery. 

The two-component analysis quantifies the ability of the model 
to predict both the correct slip directions along faults and to gen- 
erate the azimuth distribution of the faults. However, we would 
like to specifically investigate the potential of shell tectonics to 
produce the observed azimuth distribution. To do this, we include 
only those faults with slip directions that are accurately predicted 
at their current locations. We compute the probabilities of forming 
the observed azimuth distributions in each of the latitude/longi- 
tude bins from Fig. 1, using the offset accumulation rates predicted 
by shell tectonics. To compute a reference probability, we again 
use a model that assumes exclusively randomly-oriented faults. 
As described in Section 5.1, the regional probability is the product 
of the probabilities for each observed fault based on the offset 
accumulation curves. The random model probability is (l/360) n , 
but now n is different for each region. This approach allows us to 
identify regions in which the shell tectonics model performs espe- 
cially well or especially poorly. 

In Fig. 10, we show the regions in which faults were mapped by 
Sarid et al. (2002) with a color bar that indicates the ratio of the 
shell tectonics probability to the random model probability. A va- 
lue greater than 1 means that the observed azimuth distribution 
can be better reproduced by the shell tectonics model than by a 
random model. In this figure, the number of correctly predicted 
faults is listed in each region. There is only one region in which 


Fig. 11. (a) A complex tectonic region in which the azimuth distribution is better 
explained by a model that assumes random orientations than by the relative offset 
accumulation rates predicted by the shell tectonics model. The shell tectonics 
model predicts that faults with azimuths of 30-80° would form offsets at a reduced 
rate compared with faults at other azimuths, perhaps leading to these azimuths 
being under-represented in the observations. It is possible that these faults formed 
in a different location or from a process other than tidal shear stress. Faults with slip 
directions and azimuths that are well-explained by the shell tectonics model are 
shown in white; faults with black dots have azimuths that are much less likely to be 
observed according to the model, (b) A fault with an azimuth that is inconsistent 
with the shell tectonics model at its present location. The fault is connected to a 
dilational band and likely formed as a result of extension along the east-west 
trending portion (to the left in this image, as marked by the black arrows), rather 
than tidal shear stress. 

there is more than one fault and the observations support the ran- 
dom model over the shell tectonics model: longitude 210°W, lati- 
tude 15°S. Based on the relative offset accumulation rates in this 
region, we would not expect to observe offsets along faults with 
azimuths between 20° and 80°. However, 12 out of the 27 faults in- 
cluded in the analysis of this region have azimuths in this range. 
The region is shown in Fig. 11a. Faults with slip directions and azi- 
muths that are well-explained by the shell tectonics model are 
marked in white; faults with black dots have azimuths that are 
much less likely to be observed based on the model (mapping 
based on Sarid et al., 2002). 

This region is complex; geologic mapping and reconstructions 
indicate that both extension and compression have contributed 
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to the formation of bands in this area (Sarid et al., 2002; Patterson 
et al., 2006). It is possible, therefore, that some faults in this region 
were not formed as a direct result of tidal shear stress. For exam- 
ple, in Fig. lib, we show one of the faults in more detail. This fault 
connects with a dilational band shown in white. As the band 
opened along the east-west portion (as indicated by the arrows 
in Fig. lib), it would have led to right-lateral slip along the 
north-south portion. It appears that this fault was created as a 
by-product of extension, and our model correctly indicated that 
it was unlikely to have been produced by tidal shear stress. More 
in depth examination of this region may reveal that many of these 
faults were similarly formed. This would enhance our understand- 
ing of the complex tectonic evolution of this region and allow more 
accurate testing of tidally-driven strike-slip fault models. 

5.3. Additional implications 

Using the shell tectonics model and either value of the friction 
coefficient, we are able to generate predictions that fit the ob- 
served pattern of strike-slip faults on Europa. With obliquity 
(8 = 1°), the differences in the leading and trailing populations 
can also be explained. Thus, consistent with the findings of Rhoden 
et al. (2011), polar wander is not necessary to explain strike-slip 
fault observations when the effects of obliquity are included in cal- 
culations of the stress field. Also consistent with Rhoden et al. 
(2011), we find that the formation timescale of strike-slip must 
be fast relative to the NSR timescale in order for the differences 
in the leading and trailing hemispheres to be preserved. As shown 
in the “Sum” column of Fig. 4, rapid NSR should smear out the lon- 
gitudinal variations in slip direction leading to mixed LL and RL 
faults at most latitudes and all longitudes.Tidal stresses that in- 
clude obliquity are typically about 100 kPa. If faults are only al- 
lowed to slip when in tension (such as in the tidal walking 
model), the compressive stress generated by the overlying ice 
would exceed this value below about 100 m, and faults would be 
clamped throughout the entire orbit. In the shell tectonics model, 
as well as the SKP model, a fault can slip in both tension and com- 
pression as long as the accumulated shear stress is larger than the 
Coulomb failure threshold. Using our best-fit shell tectonics model 
(8 = 1°, (p = 220°, /= 0.2), faults could still slip with a compressive 
overburden stress corresponding to a depth of 400-500 m depend- 
ing on location. Hence, based on the assumptions in the shell 
tectonics model, we can only really assess the behavior of this 
uppermost portion of the ice shell. To expand our analysis to the 
whole ice shell, we would need to integrate the effects of overbur- 
den stress and depth-dependent rheology into calculations of the 
timing and amount of slip with depth. 

The statistical analyses presented in Sections 5.1 and 5.2 indicate 
that the shell tectonics model can better reproduce the slip direc- 
tions and azimuth distributions of the observed fault population 
on Europa than the SKP-type model. However, it is important to 
point out that we have used a simplified parameterization of the 
model presented by Smith-Konter and Pappalardo (2008). Our work 
tested the underlying assumptions of how fault slip is generated, 
which are different between the two models. Due to these simplifi- 
cations, we cannot fully assess the potential of the model presented 
by Smith-Konter and Pappalardo (2008) to produce offsets on Enc- 
eladus or Europa. We can conclude that the fault slip behavior they 
assumed is less likely to generate the observed fault population on 
Europa than the fault slip assumptions in the shell tectonics model. 

6. Conclusions 

Shell tectonics provides a mechanical framework to evaluate the 
response of faults within an elastic shell to periodic tidal stress. 
Shell tectonics explicitly accounts for the release of accumulated 


stress due to slip on the fault when computing subsequent fault 
stress and slip. Stress relaxation in the shell is also included explic- 
itly. Shell tectonics includes a more physical treatment of fault 
mechanics than the tidal walking model (Hoppa et al., 1999) and 
can operate much deeper into the shell because it uses a Coulomb 
failure criterion to determine when and if failure occurs. However, 
both the shell tectonics and tidal walking models assume that, in 
general, slip on a fault will increase as shear stress increases and 
will decrease as the shear stress decreases. The model of Smith- 
Konter and Pappalardo (2008) relies on a different assumption 
when evaluating fault slip. To test this fundamental difference in 
assumptions, we developed a simplified version of the SKP model 
and tested both the shell tectonics and SKP-type models against 
the population of observed faults on Europa from Sarid et al. (2002). 

When applied to Europa, both the tidal walking (Rhoden et al., 
2011) and shell tectonics models can produce global predictions of 
slip direction that display the same pattern as observed strike-slip 
faults. Consistent with previous work (Rhoden et al., 201 1 ), we find 
that an obliquity of order 1° can alter tidal stress such as to repro- 
duce the differences in strike-slip populations observed in the 
leading and trailing hemispheres. The shell tectonics model accu- 
rately predicts the slip directions on 75% of the observed faults 
from the survey by Sarid et al. (2002). For comparison, the SKP- 
type model can fit the slip directions of 50% of observed faults. 

Using the predictions of relative accumulation rates from the 
two models allows us to determine whether the azimuth distribu- 
tion of the observed faults can be explained in addition to their slip 
directions. This analysis shows that both models require a popula- 
tion of faults that formed at random orientations (50% for the shell 
tectonics model and 75% for the SKP-type model) in order to ex- 
plain the slip directions and azimuth distribution of the observed 
faults. The probability of forming the observed fault population is 
always higher with the shell tectonics model than the SKP-type 
model. 

We also determine the regional probabilities of forming the ob- 
served fault population in accordance with the relative accumula- 
tion rates as predicted by shell tectonics and a model in which all 
faults form at random azimuths. Comparison with the azimuths of 
observed faults strongly suggests that the differences in relative 
accumulation rates have led to preferred orientations within the 
strike-slip survey because faults with larger offsets are easier to 
identify. Of the 32 regions containing strike-slip faults, the shell 
tectonics model outperforms the random model in all regions con- 
taining more than one fault, with the exception of a complex tec- 
tonic region in which many faults may have formed through 
extension and compression. Further examination of this region as 
well as investigation of the role of overburden and other depth- 
dependencies on fault slip would likely improve our understanding 
of strike-slip fault formation on Europa. Shell tectonics may also be 
applicable to the formation of, and motion along, the Tiger Stripe 
fractures on Enceladus with implications for heating and plume 
generation, as these processes have already been linked to tidal 
stress (Nimmo et al., 2007; Smith-Konter and Pappalardo, 2008; 
Hurford et al., 2007). 
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